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ABSTRACT 

We present a systematic, analytical study of geometrically thin, optically thick accretion disc 
solutions for magnetized turbulent flows, with an a-like viscosity prescription. Under the 
only assumptions that (1) Magneto-Rotational instability (MRI) generates the turbulence that 
produces the anomalous viscosity needed for accretion to proceed, and that (2) the magnetic 
field amplified by the instability saturates due to buoyant vertical escape, we are able to self- 
consistently solve the disc structure equations including the fraction of power / that is carried 
off by vertical Poynting flux (and likely dissipated outside the optically thick disc). For low- 
viscosity discs, we obtain stable high-/ solutions at low accretion rates, when gas pressure 
dominates, and unstable, low-/, radiation pressure dominated solutions at high accretion rates. 
For high viscosity discs, instead, a new thermally and viscously stable, radiation pressure 
dominated solution is found, characterized by / 1 and appearing only above a critical 
accretion rate (of the order of few tenths of the Eddington one). We discuss the regimes of 
validity of our assumptions, and the astrophysical relevance of our solutions. We conclude that 
our newly discovered thin disc solutions, possibly accompanied by powerful, magnetically 
dominated coronae and outflows, should be seriously considered as models for black holes 
accreting at super-Eddington rates. 
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1 INTRODUCTION 

The so-called standard model for accretion discs 
^Shakura & Sunyaev 1973 1 describes an optically thick, geo- 
metrically thin solution to the accretion equations in which all our 
ignorance about the mechanisms of angular momentum transport 
and turbulent viscosity is side-stepped by making the (physically 
reasonable) Ansatz that the turbulent stress tensor scales with the 
local pressure, with a constant of proportionality that we will 
denote here ass - Such an approach proved enormously productive, 
and despite the great development of numerical models in the last 
decade, standard disc modeling still remains at the very heart of 
accretion disc phenomenology and the basic link between theory 
and observations. 

The recent evolution of numerical studies, instead, has proved 
itself fundamental to shed light on the very nature of the ass 
prescription, by elucidating the crucial role of MHD turbu- 
lence for the enhanced transport properties of accretion discs 
jBalbus & Hawley 1998^ . In particular, non-radiative flows are 
ideal objects of study by means of global three dimensional sim- 
ulations (see Balbus & Hawley 2002 and reference therein). 

On the other hand, the observed spectrum of thin accretion 
discs around many compact objects suggests that a corona is also 
present <Nandra & Pounds 1994 1. Such coronae are most likely a 
by-product of the internal disc dynamics (as in the model orig- 
inally proposed by Galeev, Rosner & Vaiana 1979) and local 3- 



D MHD simulations of stratified discs <Miller & Stome 2000t . al- 
though adopting a thermodynamics which is most likely inappro- 
priate for the problem at hand, may support this idea. 

The main aim of the present work is to incorporate, in the 
simplest possible way, our current knowledge on the properties of 
turbulent MHD flows in the standard accretion disc theory. This 
will allow us, on the one hand, to extend the predictive power of 
the standard model by self-consistently including in the theory the 
magnetic properties of the flow and the generation of a magnetic 
corona by means of vertical Poynting flux; and on the other hand, 
to consider MHD turbulent discs in a regime where global, 3-D 
numerical simulations are currently unable to make predictions, i.e. 
for highly luminous, radiatively efficient thin discs. 

Thus, we assume the accretion proceeds through a geometri- 
cally thin disc, a la Shakura & Sunyaev (1973), and we employ ver- 
tically integrated equations. Furthermore, we assume that magneto- 
rotational instability (MRI; see Balbus & Hawley 1998, and refer- 
ence therein) is the primary source of the turbulent viscosity and, 
consequently, that the turbulent magnetic stresses responsible for 
the angular momentum transport in the disc scale with magnetic 
pressure. Then, the accretion disc structure is fully described once 
the relationship between magnetic and disc pressure (given either 
by gas or radiation) is established. 

The issue of the viscous scaling can be translated into the un- 
certainty about the mechanisms by, and the level at, which the disc 
magnetic field saturates. In the following, we make the working as- 
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sumption (in a sense justified a priori by the need of explaining 
the observed presence of strong coronae in many accretion pow- 
ered systems) that magnetic field is disposed of mainly by vertical 
buoyancy I Stella & Rosner 1984 Sakimoto & Coroniti 1989 1, and 
verify a posteriori its validity. Then we will show how under this 
condition, and taking into account the property that MRI grows less 
rapidly in radiation dominated discs, the scaling of magnetic pres- 
sure can be determined. This in turn specifies the closure relation 
that fully describes a coupled accretion disc-corona system. 

Furthermore, we will show that by solving the structure equa- 
tions a new solution is found, which has no analog in the stan- 
dard Shakura-Sunyaev theory, and is relevant for highly viscous 
(and highly magnetized) flows accreting at a rate of the order of (or 
above) the Eddington one. 



2 BASIC EQUATIONS 

Numerical simulations of MHD turbulent discs clearly show that 
the non-linear outcome of the MRI is a fully three dimensional, 
anisotropic turbulence, that exhibits strongly correlated fluctua- 
tions in the azimuthal and radial components of velocity (Reynolds 
stresses) and magnetic field (Maxwell stress). In accordance with 
those results, and for the sake of simplicity, we will assume that 
the vertically averaged anomalous stress is dominated by Maxwell 
stresses, and therefore 



(1) 



where Pmag = is the magnetic pressure, ko is a factor of order 
unity that depends on the relative intensity of radial and azimuthal 
field components^ . 

This in turn implies that the local heating rate in a thin Keple- 
rian disc is given by 



2 ^OCs-fmag; 



(2) 



where Qk Keplerian angular velocity, Cs = ^JPuaJp is the 
isothermal sound speed. 

The growth rate of MRI is influenced by the ratio of the gas 
to magnetic pressure, as demonstrated both by analytical studies of 
the linear regime of the instability iBlaes & Socrates 2001 1 and by 
numerical MHD simulations of radiation pressure dominated discs 
jTurner, Stone & Sano 200 2 1: both found that, due to the substan- 
tial compressibility of MHD turbulence in radiation pressure dom- 
inated discs, 

a-^K-^, (3) 



where Cg = ^ P^-^sj p is the gas sound speed. The magnetic field, 
with initial (highly subthermal) amplitude Bq, rapidly grows to 
B ~ Bqc'^''', where H is the total disc scaleheight and 



(4) 



is the mean buoyant rise-time. In order to estimate the satura- 
tion field, an estimate for the upward drift velocity hd caused 



^ It is worth stressing that such an expression for the anoma- 
lous stress can be recovered exactly from the full MHD equations 
of a turbulent (thin) disc after a cai'eful averaging procedure is 
canied on i Balbus, Gammie & Hawley 1994_ .Balbus & Papaloizou 1999] 
IBlaes 2Q02i " ' 



by buoyancy of magnetic field is needed. In the case of uni- 
form discs permeated by a net mean field, where the magnetic 
field can be modeled as a collection of flux tubes that retain 
their individuality over time scales much longer than dynami- 
cal time, it has been shown j^tella c& Rosner 1984 ( that ud — 
?;A(cg/c,)(7r/2CD)'/^(a/H)i/^Herei;A = |B|(47rp)-^''2 is the 
Alfven speed, a is a flux tube cross section and Cd a drag coeffi- 
cient, whose numerical value can be estimated to be ~ 1-10. If 
the tube cross section a were independent on the ratio of gas to 
radiation pressure, such an expression for the upward drift veloc- 
ity would lead to a final scaling of magnetic pressure with the gas 
pressure only, as already argued in Stella & Rosner (1984). How- 
ever, the applicability of a flux tube analysis to the case of MRI 
turbulent flows is not justified. The work of Stone et al. (1996) 
clearly demonstrates that magnetic concentrations, although always 
present, do not persist as coherent structures. Moreover, as already 
mentioned, the turbulence in radiation pressure dominated discs is 
compressible. Then, in the vertical fluid momentum equation any 
vertical gradient of the magnetic pressure will not be balanced by 
either radiation (because radiation is diffusive) nor gas pressure 
(because azimuthal field energy density may exceed the true gas 
pressure, see below). Vertical motions are excited, with approxi- 
mately (see e.g. Eq. 4.19 in Blaes 2002) 
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which would imply (for a mainly azimuthal field) a direct propor- 
tionality between the vertical drift velocity and the Alfven speed. 
We therefore assume 



Vu = bVA, 



(6) 



with b constant of the order of unity. 

If we then put together Eqs. J3}, J4j and J6}, we have that 

-Pmag ln(i3/Bo) OC ^Ptot-Pgas, whcrC Ptot = Pgas + Prad is 

the sum of gas plus radiation pressure. Neglecting the logarith- 
mic dependence on the initial field, this in turn suggests that, in 
MRI dominated turbulent flows, the field will saturate at an ampli- 
tude such that ( Taam & Lin 1984 , Burm 1985 , Szuszkiewicz 1990l 
lMerlonr& Fabian 2002at : 



mag 



■ ao 



where we have introduced the constant 

1/2 

1 / -rtot \ 

ao 



- gas 



(7) 
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where /3 = Ptot /Pmag is the usual plasma parameter. 

It is worth noting here the relationship between the constant 



ao and the Shakura-Sunyaev coefficient ass 
is not constant in our theory). We have 



ao 



^ ass 



/Ptot (which 



(9) 



The condition that the turbulence be subsonic (va ^ Cs) implies 
ass ^1/2 (or, equivalently, /3 > 2), and from ^9) we see that ao 
can still be larger than unity, when radiation pressure dominates in 
the disc. 

The magnetic flux escaping in the vertical direction may dis- 
sipate a substantial fraction of the gravitational binding energy 
of the accreting gas outside the optically thick disc, with obvi- 
ous deep implications for the spectrum of the emerging radiation 
jHaardt & Maraschi 1991t . How much of this flux really escapes 
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into the low-density environment above and below the disc (in the 
so-called corona), depends on the details of vertical transport of 
magnetic flux tubes and vertical stratification of the disc, which 
is ignored in our model. Nonetheless, we can safely assume that 
the fraction of power dissipated in the corona is determined by 
the ratio / of the vertical Poynting flux to the local heating rate 
Q+ ISvensson & Zdziarski 199^. The simplest way to estimate 
the vertical Poynting flux is to just assume that 



together with 



Fp ~ HD-Pmag. 

This translates into (see Eq.|2j 



+ 



-Prad 
-Pgas (-R) 



-1/4 



(10) 



(11) 



In the following, we will assume that the numerical factor (of 
the order of unity) fci = 3fco/26 = 1. Differences in its exact value 
will affect the numerical value of the critical viscosity parameter 
ao (see below), but won't affect the nature of the solutions we find. 

Thanks to the closure relation <11> . we are now able to build 
a self-consistent accretion disc-corona solution. In order to do so, 
exactly as in the case of the standard Shakura-Sunyaev solution, 
we have to solve simultaneously four conservation equations (con- 
servation of vertical momentum, mass, angular momentum and 
energy) and supplement them with the equations of state for the 
pressure and opacity as functions of density and temperature. The 
only modification with respect to the standard derivation is the 
reduction of the locally dissipated energy by the factor (1 — /) 
iSvensson c& Zdziarski 1994 1, with / given by Eq. <1 1> . and the 
adoption of the viscosity law 

For simplicity, we present here the solutions in Newtonian ap- 
proximation, as in Shakura & Sunyaev (1973), with an efficiency 
eo = 0.082 appropriate for non-rotating black holes (a full gen- 
eral relativistic analog of the equations below can be found in 
the appendix of Merloni Sl Fabian 2002b). We concentrate only 
on the innermost regions of the disc, where opacity is dominated 
by electron scattering (kcs ~ 0.4 cm'^ g~^)- Introducing dimen- 
sionless units for the central mass (m = M/Mq), accretion rate 
(m = eoA'/c^/I/Edd) and distance from the center (r = R/Rs = 
Rc^/2GM), we obtain, for the radiation pressure dominated so- 
lution, the following expressions for disc scaleheight, pressure (in 
dyne), density (in g cm both calculated at the disc mid-plane, 
and temperature (in K): 



h 



H 



llmJ(r)(l-/) 

-8/9^-8/3, 



^4/5 



Rs 

P ~ 4.6 X 10^«(aom)-«/V-«/^(mJ(r))«/^(l - 
p ~ 1.9 X 10-*(aom)-«/V^/^(mJ(r))-^''/'-'(l - f)-^ 
T ~ 2.1 X lQ''{aom)-^''\-^'^{mJ{r)f''' (12) 
supplemented with the closure equation for /: 
4ag - f 



= 7.3 X 10^(Qom)'''V-^''3(mJ(r)) 



N16/9 



(13) 



Analogously, in the gas pressure dominated part of the disc, 
we have: 

h ~ 2.3 X 10-2(aom)-i/iV2i/2"(mJ(r))^/5(l - f)'^''' 
P ~ 2.7x 10^«(aom)'«/"r'^^/^°(mJ(r))^/^(l-/)-^/" 

p ~ 19(aom)-^/iV-3^/2''(mJ(r))2/5(l-/)-3/i° 
T ~ 8.0 X 10**(aom)-'/V-^/'''(mJ(r-))^/^(l-/)'/'* (14) 



4ag- 



= 4.7 X lQ\aomf^^\-^^'^\mJ{r)f'\ (15) 
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In the above formulae, the function J(r) = 1 — -\/ Vin/r, 
with Tin — 3 corresponding to the innermost stable circular orbit 
of a free particle around a non-rotating black hole, describe the no- 
torque at the inner boundary condition. 



3 STUDY OF THE SOLUTIONS 

To find all the family of possible solutions, we proceed in the fol- 
lowing way: for each set of the parameters (m, m, r and Qq) we 
find the roots of both Eq. <13> and <15> in the interval < / < 1. 
We then substitute the obtained values of / into the appropriate 
expressions for the gas and radiation pressure and finally discard 
all the non-consistent solution (i.e. those found with the formula 
for gas pressure dominated discs that gives P^ad > -Pgas and vice 
versa). 

The value of the viscosity parameter ao = fci/2 = 0.5 
is a critical point for the set of equations describing the coupled 
disc-corona system, in that it separates two qualitatively differ- 
ent regimes. In the next two sections we will examine these two 
regimes, beginning with the low viscosity case. 



3.1 Qo < 1/2: an almost standard solution 

It turns out that for ao < 1/2 each set of parameters de- 
fine a unique solution. Its properties have already been discussed 
( Merloni & Fabian 2002a i, but we present them here in a more gen- 
eral framework. The fraction of disc power dissipated in the corona, 
/, tends to its maximum value, V2ao, when gas pressure domi- 
nates (low accretion rates), and decreases as the accretion rate in- 
creases and radiation pressure becomes more and more important. 
This is shown in Fig.Q where, for m = 10 and ao = 0.25, we plot 
the radial dependence of the coronal fraction for different values of 
the accretion rate. 

As in the standard case, radiation pressure dominated part of 
the disc are thermally and viscously unstable, although the com- 
bined effects of the modified viscosity law and the stabilizing 
effect of the corona, reduce the extent of the unstable region of the 
parameter space jSzuszkiewicz 199'0l . 



3.2 Qo > 1/2: a new stable solution at high accretion rates 

The nature of the solutions changes qualitatively if ao > 1/2. Be- 
cause in gas pressure dominated parts of the disc / ^ \/2ao, there 
can be no solutions there with / < 1. That is to say that the con- 
dition ao < 1/2 should hold in gas pressure dominated discs, and 
these high viscosity solutions can only exist if radiation pressure 
dominates. 

For each value of ao > 1/2 the accretion rate m(r, /) 
(Eq. Us) has a minimum for / £ (0, 1). Thus, solutions exists 
only for accretion rates larger that a critical value, rhcr{r) = 
min/g(o,i) [m(/, r)] = m(r,/cr). This corresponds to a critical 
value, /„■, of the fraction of power vertically transported by Poynt- 
ing flux. For m(r) > rhcrir) two radiation pressure dominated so- 
lutions are found for each set of parameters, one with / < fcr and 
the other with / > /„ (see Fig.|2}. The first is a "standard" (un- 
stable) one, characterized by a low value of the fraction of power 
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Figure 1. Radial profiles of the fraction of power transported vertically by 
Poynting flux, / for m = 10 and ao = 0.25 and different values of the ac- 
cretion rate: m = 0.003, 0.01, 0.03, 0.1, 0.3, 1, 3. The disc structure has 
been computed using Newtonian approximation for non-rotating black hole. 
Filled triangles represent gas pressure dominated solutions; crosses radia- 
tion pressure dominated ones, while empty circles mark the solutions where 
cooling by advection should be taken into account (slim disc solutions). 
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Figure 2. Radial profiles of the fraction of power transported vertically 
by Poynting flux, / for m = 10 and ag = 0.75 and different values 
of the accretion rate: rh = 0.2, 0.4, 0.8, 1.6, 3.2, 6.4. The disc structure 
has been coinputed using Newtonian approximation for non-rotating black 
hole. Filled squares represent stable, radiation pressure dominated solu- 
tions; crosses unstable ones, while empty circles mark the solutions where 
cooling by advection should be taken into account (slim disc solutions). The 
dotted horizontal line mark the critical value fa ~ 0.689 that separates the 
two solutions. 



transported vertically by Poynting flux, and is just the continua- 
tion of the radiation pressure dominated branch discussed in the 
previous section at high values of uq. The second one, which is 
discussed here for the first time, appears due to the feedback effect 
of the closure relation JUL It is a (mildly) radiation pressure dom- 
inated solution endowed with a powerful corona, in that / ^ 1 as 
the accretion rate increases. 

The values of mcr(r) and can be found by solving 
drn{r, f) /df = 0, or, equivalently, the equation f"' — 12ao/ + 



Sqq ~ 0. For every r, re obtain /„■ ~ 0.7417 for ao = 0.5, 
quickly converging towards fcr = 2/3 for increasing qq. The 
critical accretion rate, though, depends on r. The lowest possi- 
ble value of the critical accretion rate, that can be simply denoted 
rhcr, appears at r ~ 4.64 (which is the radius where the function 
j('^^i6/9j,-7/3 jj^^ ^ minimum). We have numerically computed 
rhcr as a function of qq. The relation can be approximated by the 
analytic expression 

-1/8 



rhcr ~ 0.34aom 



(16) 



which is accurate to better than 3% for qq > 1 and to better than 
10% for ao > 0.6. The above formula also makes manifest the 
main differences that would be found if we chose supermassive 
black holes instead that stellar mass ones. Namely, the critical ac- 
cretion rate separating gas from radiation pressure dominated solu- 
tions and that determining the existence of our new high viscosity, 
high-/ solution would be reduced by a factor m^^''*. The qualita- 
tive properties of the solutions would be left unchanged. 

Fig-IHshows, for m = 10 and ao = 0.75, the radial depen- 
dence of the coronal fraction for different values of the accretion 
rate. The two solutions appear for m>0.19 and show an oppo- 
site trend as the accretion rate increases: while in the first case / 
decreases and the disc becomes less dense and more radiation pres- 
sure dominated, in the second case / increases, the disc becomes 
denser and less radiation pressure dominated. 

It is worth studying the stability property of this newly discov- 
ered solution. Thermal stability requires iPiran 197 8 J 

'd\nQ+\ /91nQ. 



/ ainQ+ \ / ginQ_ \ 
V dlnT )r^\ dlnT ) ., 



(17) 



Wehave(ainQ_/91nT) = 4 and, from GI}. <3+ « m(l-/) oc 
7-9/2(1 _ J). Therefore, the stability criterion becomes: 



/ ainQ+ \ _ /91nO_\ ^ 9 / 
V ainT Jr \ dlnT Jr^ 2\ 



1 + 



dlnT J ^ 2 

4/(4ag - f) 



1 

2 + (7 



d{l-f) 
dm 



< 0. 



(18) 



Ufal + Sag) 

Because the numerator of the second term is always positive for 
ao > 1 /2, we have that for / > fcr (which is a pole for the second 
term in the above equation) indeed the criterion is satisfied, while 
the opposite is true for / < /„ • On the other hand, viscous stability 
requires that {dM/dT.) > 0. We have E oc pff oc m"^/'-*(l - 
/)~^, and the stability condition can be written as 

4/(4ag 



2 / 1 



n 

12fal + Sag 



< 0, 



(19) 



which is satisfied when thermal stability condition is also satisfied. 

The very important fact that our new discovered high-/ so- 
lutions are both viscously and thermally stable may be an indica- 
tion that, whenever allowed to do so, the system will chose those 
solutions over the low-/ ones. The real physical state of a high 
viscosity disc at high accretion rates, though, is likely more com- 
plicated and inherently time dependent. Indeed, the very fact that 
for ao > 0.5 no gas pressure dominated solutions exist, make 
us believe that in the outer part of the disc magnetic turbulence 
will self-regulate to keep the viscosity parameter below its criti- 
cal value. Thus, the value of the viscosity parameter may not be 
constant throughout the disc, but more likely increases inwards, as 
also suggested by high-resolution numerical MHD simulations of 
the inner parts of non radiative discs ^Reynolds & Armitage 2001] 
Krolik & Hawley 2002 (. For any given radial profile of the viscos- 
ity parameter ao, the radius at which it reaches its critical value will 
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therefore play an important role in determining the physical prop- 
erties of the disc. We may speculate that in a full time dependent 
solution both the viscosity parameter and the local accretion rate 
will adjust in order to avoid discontinuities in the flow there, but 
the exact nature of the final disc structure is impossible to predict 
from our stationary solutions. 



4 DISCUSSION 

Beside the thin disc approximation, which is straightforward to ver- 
ify a posteriori for all our solutions but the low-/, high-rh ones (as 
in the standard SS case), our results are rooted in two further as- 
sumptions. The first is that the growth rate of MRI is given by J3j 
in radiation pressure dominated discs, the second is that buoyancy 
is the main mechanism of field saturation in turbulent magnetized 
discs. Let us examine these two issues one at a time. 

As discussed in Blaes & Socrates (2001), the applicability of 
Eq. J3} depends on the capability of radiative diffusion to destroy 
radiation sound waves on the length scale of the MRI, and requires 
that the azimuthal magnetic field component dominates the vertical 
one, as should be the case in thin discs with Keplerian velocity pro- 
files (see also the numerical simulations of initially random fields in 
Hawley, Gammie & Balbus 1996). Let us define the Alfven speed 
associated to the vertical field component, uaz = BziA-Kp)'^^'^ , 
and the (small) parameter 5 = vaz/va- The diffusion wavenumber 
for Keplerian discs is (Bl aes & Socrates 2001j : 



fcdiff 



1/2 



(20) 



Then, Eq. should be valid if the characteristic wavenumber 
of the instability, Qk/vaz is of the order of (or larger than) the 
diffusion wavenumber fcdift . By using the expression fl2\ for the 
density, <1 H . and the definition of vaz, we obtain the following 
condition for the validity of our assumption: 



kdisfHS ~ 65/m 



Noting that max [J(r 



1/18 



mJir) 
Qo 



4/9 



-21/36^ <^ 



,4/9^-21/361 



: 0.2, this translates into 



1.5/mj 



1/18 



(-) 



4/9 



<1, 



(21) 



(22) 



which is always satisfied for the low-/ radiation pressure domi- 
nated solutions, while, for the / ~ 1 solutions, it breaks down 
for very high accretion rates m ^ 0.55ao?n~^^**(5Zi''^ — 
L6mcr<5Zi''*- It is clear that in order for our viscosity scaling Q 
to be applicable to the high-/ solutions at super-Eddington rates, 
higher azimuthal field strengths (and low 5) are required. However, 
if m grows too large, and Eq. <22t is not satisfied, then the magnetic 
stresses, and the magnetic pressure, will be proportional to the total 
pressure. The condition of subsonic turbulence would then imply 
again ao < 1/2. The fraction of power transported vertically by 
Poynting flux, / will tend again to its maximum value ~ \/2cto, as 
in the gas-pressure dominated solutions. Therefore, it seems plausi- 
ble to envisage that the high-/ solutions could extend to arbitrarily 
large accretion rates, regardless to the viscosity law and the de- 
tailed properties of the MRI field, as long as qq is sufficiently high 
in these regimes. 

The second main assumption that was used in the derivation 



We adopt here the notation A = IQ^ Ax- 



of our solutions was to consider buoyancy the main mechanism 
by which the disc gets rid of the magnetic field. In fact, two other 
competing mechanisms should be considered: turbulent diffusion 
and dissipation. 

Let us first compare the typical buoyancy timescale {tb — 
H/2bvA) with the turbulent diffusion one (td — /v). The kine- 
matic viscosity coefficient u can be recovered from the standard 
viscous form of the stress to obtain v = 2tr4, / {?>pQ.-k)- We then 
have 



66cs 6 



(23) 



Therefore, for any physical self-consistent, / < 1 solution, buoy- 
ancy is always a more efficient mechanism of field escape than field 
diffusion, unless fci = 3fco /2fe is substantially larger than unity.. 

On the other hand, the field may dissipate by cascading down 
the subsonic turbulence and leave its energy inside the disc. The 
typical timescale for this process is set by the coherence time of 
the turbulent flow, that we have assumed to be Tt = 1/a ~ 
VA/{cgi^K)- This would give 



Tt 



46an 



(24) 



which implies that our working assumption is justified only for 
large viscosity solutions. In practice, if ao is small, the actual max- 
imal fraction of power dissipated into the corona will be further 
reduced by a factor 1 + 0.25/feao- However, the high viscosity, 
high-/ solutions discussed in section 3.2 will be little affected by 
dissipational losses of energy inside the disc. 



5 CONCLUSIONS 

We have presented a full analysis of thin accretion disc solutions 
under the assumptions that the magnetic field generated and ampli- 
fied by MRI (which is at the origin of MHD turbulent stress in the 
disc) saturates mainly due to buoyant transport in the vertical direc- 
tion. From these assumptions, a viscosity prescription is derived, 
in which magnetic pressure, and therefore turbulent stresses, scale 
proportionally to the geometric mean of gas and total (gas plus ra- 
diation) pressure. This prescription, in turn, allows us to uniquely 
determine, for every set of the parameters (m, m, r and ao) the 
fraction of the accretion power / which is transported vertically by 
Poynting flux. For low values of the viscosity parameter the solu- 
tion is unique, with / tending to its maximum value, \/2ao, when 
gas pressure dominates (low accretion rates), and decreasing as the 
accretion rate increases and radiation pressure becomes more and 
more important. For 2ao > 1 there are no solutions for gas pres- 
sure dominated regions, while in the radiation dominated part of 
the flow two solutions are possible for every value of the accretion 
rate. The first is a "standard" (unstable) one, with / ^ 1, while 
the second has / ^ 1 (corona dominated) as the accretion rate 
increases and is thermally and viscously stable. 

Two main results from the present study are important for 
the interpretation of observed spectral energy distributions of 
accretion-powered systems. The first is the property of low-ao so- 
lutions to be more corona-dominated (i.e. to have the largest pos- 
sible value of /) at low accretion rates. This has been discussed in 
Merloni & Fabian (2002a) in the context of low-luminosity black 
hole. It is indeed an observational fact that the relative strength of 
the hard, optically thin, power-law-like component in their X-ray 
spectra increases with respect to the optically thick, quasi-thermal 
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one as rh decreases. The second is the existence of a new high- 
/ solution at high accretion rates for high values of the viscosity 
parameter. 

As in the case of corona-dominated solutions at low m, it 
is plausible to speculate that the large Poynting fluxes in the ver- 
tical direction that characterize these solutions are to be associ- 
ated with the generation of powerful coronae and, possibly, pow- 
erful outflows/jets, too t Mer loni & Fabian 2002a> . Indeed, there 
are observational hints that black holes believed to accrete at 
Eddington or super-Eddington rates bear many similarities with 
their low-rate counterparts. This class of objects may include 
galactic black holes in the Very High state (see Done 2002, 
and references therein), whose X-ray spectra show a power-law 
component, likely of non-thermal origin i Pout anen & Coppi 1998| 
[Gierliriski et al. 1999 1, often as strong as the quasi-thermal one; 
broad line radio galaxies, where the weakness of X-ray reflec- 
tion features has been interpreted as evidence for highly ion- 
ized discs surrounded by powerful coronae at high accretion rates 
jBallanty ne, Ross & Fabian 2002'>; powerful blazars, where it has 
been shown I Ghisellini & Celotti 2002 ( how the relativistic jet in 
these sources may dominate the total output power at very high 
(inferred) m. Furthermore, some Ultra Luminous X-ray sources 
(ULX; Makishima et al. 2000) have shown power-law like spec- 
tra, analogous to those of black hole candidates in their very high 
state, at inferred accretion rates of a few times the Eddington rate 
for a 30Mq black hole jKubota, Done & Makishima 2002t . 

In order to assess the relevance of our new solution in all the 
above cases, vertically integrated models should be abandoned, and 
the complicated problem of determining the disc structure in the 
vertical direction, together with the disc-corona boundary, tack- 
led. This would allow us to establish the relationships between the 
(vertically integrated) parameter /, and the observationally deter- 
mined fraction of power dissipated in the corona. In particular, it 
is only by carefully studying the disc-corona transition region that 
will be possible to assess how super-Eddington (if at all) our new 
high Qo solutions can be. In principle, if a high fraction of ac- 
cretion power could be dissipated in the corona without disrupt- 
ing the optically thick disc, there would be no limitation to the 
total luminosity that the source could emit, because the radiation 
is shining on the disc from outside. Such a possibility to obtain 
super-Eddington luminosities is physically different from the one 
recently proposed by Begelman (2001; 2002), although both de- 
pend on the existence of magnetized radiation pressure dominated 
thin discs. In fact, the non-linear analysis of Gammie (1998) and 
Blaes & Socrates' (2001) photon bubble instability carried on in 
Begelman (2001) is applicable to the low-/ radiation dominated 
solutions we have presented here, which is much more strongly ra- 
diation pressure dominated than the high-/ ones, while the thinner 
discs of the high-/ solutions have almost equal radiation and gas 
pressures, and therefore should not be too strongly inhomogenous. 
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Finally, we point out that the existence of multiple solutions 
for the same set of parameters (high qq case), may give rise to in- 
teresting hysteresis effects during spectral transition caused by long 
term variation of the accretion rate (Maccarone & Coppi 2003^ . A 
full time-dependent analysis of our solution is needed in order to 
make definite predictions in this sense. 



